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Abstract 

A multispecies artificial ecosystem is formulated 
using cellular automata with species interactions 
and food chain hierarchy. The constructed finite 
state automaton can simulate the complexity and 
self-organized characteristics of the evolving multi- 
species living ecosystems. Numerical experiments 
show that a small perturbation or extinction event 
may affect many other species in the ecosystem in 
an avalanche manner. Both the avalanches and 
the extinction arising from these changes follow a 
power law, reflecting that the multispecies living 
ecosytems have the characteristics of self-organized 
criticality. 

Citation detail: X. S. Yang, Characterization 
of multispecies living ecosystems with cellular au- 
tomata, in: Artificial Life VIII (Eds. Standish, 
Abbass, Bedau), MIT Press, pp. 138-141 (2002). 



1 Introduction 

The modelling of artificial life is very important 
to the understanding of the origin of life and its 
evolution, thus the studies using computer simula- 
tions are of both scientific importance and practi- 
cal applications. Since Langton (1986) relauched 
the field of artificial life, the possibly universal as- 
pects of living systems have been investigated by 
exploring the artificial chemistry acting on artificial 
molecules in terms of rule-based cellular automata 
(CA). Langton (1986) carried out numerous CA 
simulations with a parameterization scheme allow- 
ing the relationship of Wolfram's classes cataloged 
the rules that generates different classes of dynam- 
ics. Since then there have been many extensive 
studies about artificial life. However, the under- 
standing and studies of artificial life is still at very 
early stage. 

Artificial life simulations using cellular automata 
have received much of interests in the community 



of artificial life and evolutionary computing. Con- 
way's game of life has popularized this area of re- 
search. However, most of the existing simulations 
of CA is about the evolutionary and classification 
of one single species, and this has demonstrated 
the complexity and richness of simple rule-based 
systems. Interestingly, multispecies models have 
begun to attract much attention (Sole and Man- 
rubia, 1996; Amaral and Meyer, 1999). In order 
to study the ecological effect and the interaction 
among different species, a multispecies system of 
artificial life should be explored in detail. It can be 
expected that the behavior and characteristics of 
multispecies living system could be very different 
from that of single species cellular automata. 

The aim of this paper is to present some new 
results using a multspecies system of artificial life 
with cellular automata. By numerical experiments, 
we will find how the interactions among differ- 
ent species affect the evolution behavior and self- 
organization as well as species extinction events. 
This paper is orgainised as follows. In the next sec- 
tions, we describe how to construct the cellular au- 
tomata for the multispecies living systems, we then 
implement the systems and give some numerical 
simulations. Based on the numerical simulations, 
the complexity and entropy of the evolving living 
systems are measured, and the self-organized crit- 
icality is tested. The possible reason of extinction 
behavior of the living systems is also explored. The 
implications of the multispecies interactions will be 
also discussed. 



2 Rule-Based Cellular Au- 
tomata for Multi-Species 
Living Systems 

Cellular Automata (CA) is a finite-state machine 
on a regular lattice. The input to the machine is 
the states of all machines in its neighborhood, the 
change of its state is based on the rules or transition 
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functions. The states of all machines in the lattice 
are updated synchronously in discrete time steps. 
Each cell evolves according to the same rules which 
depend only on the state of the cell and a finite 
number of neighboring cells, and the neighborhood 
relation is local and uniform. One classic example 
is Conway's game of life where one takes a neigh- 
borhood consisting of the nearest 8 cells to a cell 
on a 2-D cellular automata. Its transition function 
for the local automaton for two states (1 and 0, or 
alive and not alive) is as follows: If 2 or 3 neighbors 
of a cell are alive (or 1) and it is alive at present, 
then it is alive at next state; If 3 neighbors of a cell 
are alive and it is currently not alive, its next state 
is alive; the next state is not alive for all the other 
cases. Even with these simple rules, an universal 
computational machine can be constructed on an 
infinite 2-D grid, which is capable of emulating the 
computing power of any Turing machine or digital 
computer. 

The present model is an extension of the com- 
bined version of the Sole-Manrubia model (1996) 
and Amaral-Meyer food chain model (1999). The 
next state for a cell is determined by the transition 
function in the general form 

Ci{t + l) = H[J2GijCi{t)-9], (1) 

i 

where H is the heaviside step function, and 6 is 
the threshold. Gij of the infiuence of cell j on cell 
i, Gij > or < if j is the food or predator for 
cell i. The summation is over the nearest neighbor- 
hood. The multispecies cellular automata works in 
the following way: 

1. The CA consists of M interacting species, 
when M = 1 it degenerates into the classic 
Conway's game of life and obeys the same sim- 
ple rules. 

2. Different sepcies interacts with each other in 
a way of food chain hierarchy. Each species is 
label as a hierarchy level and for simplicity, we 
use the level i for species i with species i — 1 as 
the pray or food for species i and species i + 1 
as its predator. 

3. In the nearest 8 cells of to a cell of species i, 

if the number of predators is greater than the 
number of prays, and the cell is currently alive, 
then its next state is alive; If the number of the 
prays is more than the number of predators, 
then its next state is alive whether its present 
state is alive or not; If the predators are in the 
same number as the pray, the present state 
does not change. 



4. If there is no predator or pray in the 8 neigh- 
borhood cells, then the transition function is 
the same as the single species Conway's cellu- 
lar automata. 

In the numerical simulations, the 2-D lattice is 

randomly initialized with the highest population 
of the lowest level of species i = 1. For a regu- 
lar 100 X 100 lattice with N=256 species initial- 
ized this way, wc can investigate the complexity, 
self-organization, species interaction and avalanche 
events. 

To simulate the effect of species interaction and 
extinction, a small probability of extinction rate is 
also introduced into some species, the infiuence of 
the extinction of one species on the other species in 
the food chain model can be studied in some detail. 
In addition, by introducing some extra population 
in one species, one can control some other species in 
the food chain model and give important implica- 
tions on the ecological eff'ect in reality such as the 
rat control. A number of numerical experiments 
have also been investigated. Some trend and self- 
organized criticality is analyzed from the simulated 
results. 

3 Numerical Experiments 

By using the cellular automata described in the 

above section, wc can simulate the behavior of 
multspecies living systems with multispecies inter- 
actions in a food chain hierarchy. In the rest of 
this paper, we present some results from a huge 
number of the numerical simulations and paramet- 
ric studies. A typical initial configuration is shown 
in Figure 1. 

3.1 Variation in Complexity and En- 
tropy 

For a population size N = 100 x 100 with M = 
256, the complexity of the cellular automata can 
be measured by its entropy S 

S = -Y.Pi'^ogpi, (2) 

i 

where pj is the probability of find the species i 
in the total population (Adami,1998). For a fi- 
nite population size, pi can be approximated by 
the fraction of samples rii 
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Figure 1: A random configuration of multispecies 
cellular automata 
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Figure 2: Variations of complexity versus time 
steps 

so that 

5" = -^Pilogp.. (4) 

i 

The variation of complexity of the 256-species cel- 
lular automata is shown in Figure 2. 

It is clearly seen that the complexity varies sig- 
nificantly at the early stage of the evolution pro- 
cess, then it gradually relaxes to the equilibrium at 
long time, indicating that the living systems is in a 
quasi-steady state among different species. 

3.2 Self-Organized Criticality 

Based on the pioneering work by Adami (1995) on 
the self-organized criticality in living systems us- 
ing the Tierra experiments and the subsequent 
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Figure 3: Self-organized criticality in multispecies 
living systems 

work by Newman et al (1997), there is a good 
reason to believe that self-organized criticality ex- 
ists in the artificial living systems. In order to 
test whether the self-organized criticality exists in 
the multispecies living ecosystems, we introduce a 
small perturbation to ecosystem and measure the 
population changes affected by the perturbation or 
avalanche in population and species, this is because 
the ecosystem can form between all the competing 
and yet interconnected species. The change in re- 
sponse to any perturbation will enable this living 
to relax to a new equilibrium or the self-organized 
state. By tracing the 15000 avalanches obtained 
from numerical simulations on 100 x 100 lattice 
grid, the results are shown in the log-log plot in 
Figure 3. 

This figure clearly show that the avalanches size 
distribution follows a power-law with the exponent 
7 = 2.26 ± 0.04. This implies that the muhi- 
species living ecosystem has the characteristics of 
self-organized criticality. This is consistent with the 
earlier work by Adami (1995). Our present work 
implies that the self-organized criticality could be 
universal in living systems. 

3.3 Extinction 

In the above numerical simulations, the changes 
or redistribution of populations of different species 
have been studied. There may not necessarily in- 
volve the extinction of some species. In some case, 
for example, when the initial population size for 
some species is very small and its predators are 
strong, the whole species could become extinct, and 
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Figure 4: Extinction intensities at different time 
stage 

the extinction of one species may influence other 
species in the food chain. To simulate this effect, 
we introduce an extinction probability pe = 0.01 
for a species at some level, say i = 16, in the food 
chain hierarchy. The extinction obtained from the 
numerical simulations is shown in Figure 4. 

The extinction intensities at different time stages 
(100 time steps as a stage) follow a power- law with 
an exponent 7 — 2.05 ± 0.1, which implies that 
extinction is a self-organized phenomenon. This is 
consistent with the fitting data j — 1.7 ± 0.3 from 
the fossil record (Newman and Eble, 1999), and the 
extinction is a result of external change and species 
competition. 

4 Discussion 

The multispecies ecosystems of artificial life have 
been formulated in terms of 2-D cellular automata 
of the combined Sole-Manrubia type together with 
the Amaral-Meyer food chain model. By using the 
proper interactions between different species in the 
food chain hierarchy and the transition function, 
the finite state automaton can simulate the com- 
plexity of the evolving multispecies living ecosys- 
tem. Numerical experiments show that the com- 
plexity measured by the system entropy fluctu- 
ates at early stage of evolution, then it evolves 
to the self-organized or dynamic equilibrium sate 
under certain conditions. By introducing the ex- 
tra fraction of population or a small perturbation, 
the species population may be affected at different 
scales. The avalanches arising from these changes 
obey a power-law with the exponent 7 = 2.26, re- 



flecting that the multispecies living ecosystem has 
the characteristics of self-organized criticality. 

One the other hand, in the case of one species in 
the margin of extinction, it can affect many species 
in the multispecies living systems and can lead 
some kinds of mass extinction. Numerical simu- 
lations also show that the extinction intensities fol- 
low a power-law form with the exponent 7 = 2.05, 
thus indicating the species interaction and competi- 
tion is the mechanism for the extinction under the 
action of external changes. Although the present 
modelling provides a lot of features about the com- 
plexity of the artificial living systems, much more 
work is clearly needed to investigate the system be- 
havior such as how the different extinction mech- 
anisms and transition rules affect the evolution of 
the multispecies living ecosystems. 
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